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ABSTRACT 



Aims. The aim of this paper is to study the characteristics of the stellar populations and the metallicity distribution in the Galactic 
bulge. We study the entire stellar population, but also retrieve information using only the red clump stars. 

Methods. To study the characteristics of the stellar populations and the metallicity distribution in the Galactic bulge, we compared the 
output of the galaxy model TRILEGAL, which implements the Binney et al. (1997) bulge model, with observations from 2MASS and 
OGLE-II. A minimisation procedure has been set up to retrieve the best fitting model with different stellar populations and metallicity 
distributions. 

Results. Using the TRILEGAL code we find that the best model resembling the characteristics of the Galactic bulge is a model with 
the distance to the Galactic centre R = 8.7+^3 kpc, the major axis ratios of the bar 1 : 77 : £ = 1 : 0.68±[j ^ : 0.31±g'|jJ, and the angle 
between the Sun-centre line and the bar tp = \y± i ^ i 1 . Using these parameters the best model is found for a burst of 8 Gyr, although 
it is almost indistinguishable from models with ages of 9 and 10 Gyr. The metallicity distribution found is consistent with metallicity 
distributions in the literature based on spectroscopic results. 



Key words. Galaxy: bulge - 
and C-M diagrams 



1. Introduction 



Galaxy: structure - Stars: distances - Galaxy: fundamental parameters - Stars: Hertzsprung Russel (HR) 



The central part of our Milky Way Galaxy (MWG) consists of 
the Galactic bulge (GB), the central part of the halo, and the 
Galactic centre (GC). Over the years, several studies have tried 
to understand the characteristics of these parts of the MWG, for 
example about the formation history, the rate of star formation, 
the relatio ns between different stellar components, and its 
kinematics dWvsdll 99*51) . 

Understanding the formation history of galactic spheroids 
(i.e. elliptical galaxies and the bulges of spiral galaxies) is 
crucial in understanding galaxy formation in general as these 
spheroids contain a large frac tion of all stellar mass in the local 
universe dZoccali et al.l 12003, and references therein, hereafter 
ZRO2003). Because of its proximity, our GB is the spheroid 
for which the most information is accessible. Issues that can 
be addressed in the GB are, e.g., the initial mass function, the 
distribution of stellar ages, and the 3D structure. For more 
distant spheroids, this sort of information is only accessible 
by very indirect means based on the study of the integrated light. 

Many different approaches exist to studying galaxy forma- 
tion. Often galaxy models are used to study the entire content 
of a galactic component such as the bulge dDwek et al.l IT995t 



Binn ey etal.1 119971; iFreudenreichl 1 19981 Bissantz & Gerhard 



2002; iRobin et all 120031; iMerrifieldl |2004|) . but studies of 



* Figure 9 is only available 
http://www.edpsciences.org 



in electronic form via 



the same kinds of stars can also reveal a wealth of infor- 
matio n. Examples of th i s are, e.g., the use of red clump 
stars (IStanek et all 1 19971 iPaczvnski & Stanekl 119981; [Od alski 
1 19981; iBabusiaux & Gilmorel 120051; iNishivama et alj ]2006: 
Lopez-Corredoira et al. 2007; Rattenbury et al. 2007), RR Lyrae 
stars (iFernlev et al l Il987b iFeastl 119971 lUdalskil Il998l) . masers 
(lReidetal.1 ll988USevensteretal.l Il999h. and Mira variables 
dWhitelockll992t [Groenewegen & Blommaertl2005l) . Our target 
is both to perform a calibration of the bulge geometry and to ex- 
plore the distribution of ages that can be inferred from these data. 

The paper is organised as follows. In Sect. [2] we give a 
resume on the models and parameters generally used to study 
the GB, the star formation rate, and the metallicity distribution 
of the GB. In Sect. [3] we give a brief introduction to the 
galaxy model used in this work. Section |4] is dedicated to the 
observational data used. Section [5] explains the method used 
and Sect. [6] discusses the results derived using this method. 
Section Q discusses the results using the same method if only 
the red clump stars are used. In Sect. |8]we summarise and come 
to the conclusions. 



2. The Galactic bulge 

At the moment, it is still not clear whether one should draw a 
distinction betwe en a triaxial central bulge and a bar component 
dMerrifieldll2004l). Sometimes t he bulge is referred to as the bar 
(see e.g. ISevenster et al.lll999l) and in other cases studies find 
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a bar within the bulge (see e.g. lLopez-Corredoira et al.l 120071; 
Cabr era-La vers et al.l |2007 | ). In this pap er we will adopt the 
formalism described bv lSevensterl d 1999b : we will use the term 
Galactic bulge for the Galactic component seen in the general 
direction of \£\ < 10°, we will use bar if we refer to the prolate 
or triaxial component of the GB. 

This section gives a brief overview of recent studies propos- 
ing different models and parameters describing the GB. Also the 
possible star formation rate and age metallicity relations will be 
described. 



2. 1 . Models describing the Galactic bulge and its bar 

Many different geometries can be assumed for the GB. Oblate 
spheroidal models can be a good start, providing a first order 
esti mate of the bulge luminosity and minor to major axis ra- 
tios dDwek et alj|1995l) . As became clear over the past decades, 
there is a triaxial part, the bar, which needs a triaxial model. 
There are three types of triaxial models proposed: Gaussian- 
type funct ions, exponential- type functions, and pow er-law func- 
tions ( see lDwek et al.ll 1995L and references therein) . [Pwek et ail 
d 1995b find that the best model to fit the COBE/DIRBE data 
for the bulge is a n exponentia l -type function with a modified 
spher oid included . iRobin et al.l d2003l) use the same parameter 
set as IDwek et al.l dl995h~ introduced to fit the bulge. They used 
a boxy Gaussia n type funct i on, bu t did not model the inner parts 
of the bulge. iBinnev et al.l d 19971) and also [Bissantz & Gerhard 
(2002) used a truncat ed power-law to mo del the bulge. In this 
paper, we will use the IBinnev et al.l (1 19971) formalism to model 
the number density of stars /bulge, which is given by: 



Cabrera-Lavers et al 



A- 



bulge 



(1 + fl/flo) 



and 



(i) 

(2) 



In these formulae, a m represents the scale length, «o the inner 
truncation length, and /o is the stellar density per kpc 3 of the 
bulge. The axis-ratios of the bar are given by: 1 : rj : 



C abrera-Lavers et al 



the red clump st ars in the GB c o mes fro m 
(120071) and iRattenburv etail d2007l) . 
(120071) claim two very different large-scale triaxial structures 
in the inner Galaxy. A first component is a long thin stellar 
bar Qb\ < 2°) with a position angel of 43.0° + 1.8°. The 
second component is a distinct triaxial bulge that extends to 
at least \b\ < 7 . 5° an d has a position angle of 12.6°± 3.2°. 
IRattenburv et al.l d2007l) study the red clump stars observed with 
OGLE in the OGLE-II phase and find that the bar's major axis 
is oriented at 24° - 27° to the Sun-Galactic centre line-of-sight. 
Another kind of stars that is considered a standard candle are 
the RR Lyrae stars. Using these stars typical values for Rp of 
~ 8.0 kpc have bee n found (see e.g. iFernlev et al.lll987l: iFeastl 
U997UUdalskilll998l) . 

Other types of stars that have bee n studied to characterise the 
GB are e.g. Mira variables (see e.g. Groenewegen & Blommaert 
(120051 R = 8.8 kpc and = 43°) and IWhitelockl dl992l 
1 : 77 : £ = 1 : 0.25 : 0.25, <p = 45°)), H 2 masers 



(see iReid et al l 1 19881 R = 7.1 kpc), and OH/IR stars (see 



ISevenster et al 



1999, a m = 2500 pc and = 44°). 



Studying one type of star is not the only method to obtain 
information on the characteristic s of the Galacti c bulge and 
distance to the Galact ic centre. lEisenhauer et al.1 (120031) and 
Eisenhaue r et al.l (120051) studied stars orbiting the central black 
hole and found values for Rg of respectively 7.94 + 0.42 kpc and 
7.62 ± 0.32 kpc. Recently, iGhez et al.1 d2008) and lGillessen et al.l 
(2009j) used additional data to find, respectively, 8.0 + 0.4 
kpc and 8.33 + 0.35 kpc. Also surveys like 2MASS and the 
surface brightness maps of COBE/DIRBE were used. Studies 
based on these databases find values for </> around 20° (see 
e.g. IDwek et al.1 119951; IBinnev et al.l 1 1 997b iBissantz & Gerhard! 
120021: lLopez-Corredoira et al.l [20051). scale lengths a m between 



1900 pc and 2800 pc (see e g. Binnev et"aTlll997t iFreudenreichl 



119981: IBissantz & Gerhard! |2002|). and an axis-ratio o f 1:0.5- 
0.6:0.4 dMerrifield|l2004t lLopez-Corredoira et alll2005l) . except 
lDweketal.1 dl995[) who find low er values for the axis-ratio 
(1:0.33:0.23). IRobin et all (120031) used HIPPARCOS data to 
model and retrieved axis-ratios of 1:0.27:0.27 and an angle <p of 
11.1° + 0.7°. 



2.2. Parameters describing the Galactic bulge 

The typical parameters needed to describe the GB (see previous 
section) are the scale length (a m ) and inner truncation length («o) 
of the bulge, the axis-ratios (1 : 77 : £) of the bar, the angle be- 
tween the Sun-centre line and the major axis of the bar (0), and 
the scaling parameter (Jo). Of course also the Sun's distance to 
the Galactic centre (Rq) is an important parameter in this context. 



The Milky Way Galaxy consists not only of stars, there 
is also an amount of gas pre sent. The locatio n of Hi gas 
in the galaxy (see Fig. 2 in Merrifield 2004) immediately 
reveals the non-axisymmetric distribution of the gas. This non- 
axisymmetric distribution can also be seen when o ne measures 
the v elocities of the interstellar CO mol ecule dDame et al.1 
120011) . Using the distribution of the Hi gas Merrifield (200' 
retrieve an angle <p of 25°. 



One of the most studied types of stars in the GB are 
the red clump stars. The red clump stars are considered as 
standard candles to determine distances. They can also be 
used to estimate the axis ratios of the bulge and the angle <p. 
Using these red clump stars, value s between 7.4 and 8.4 kpc 
for Rq have b e en found (see e.g. iPa czynski & Stanek 1998; 
UdalsB 119981; IBabusiaux & Gilmord 120051; iNishivama et al l 
20061) . Values for <p derived from red clu mp stars are s i tuated 
either around 2 0° or around 45° (see e.g. | Stanek et al.l (1 19971 
<f> = 20° - 30°), IBabusiaux & Gilmorel (|2005j = 22° t 5.5°), 
and lLopez-Corredoira et al.1 d2007l <p = 43°)). Recent work on 



Table Q] gives an overview of these studies and the parame- 
ters retrieved. Concerning the angle between the Sun-centre line 
and the bar </>, Table[T]clearly shows two groups of values: a low 
value around 20° and a higher value around 45°. According to 
ISevenster et al.l d 1999b the lower values found for cf> arise when 
the longit ude range used is too narrow or whe n low latitudes are 
excluded. iGroenewegen & Blommaert (2005) point out that this 
could also be due to the fact that these studies trace different 
populations, which may be distributed differently. 
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Reference 




a m a 1 : r/ : ( 





Based on 




[kpc] 


[kpc] [pc] 


n 





Fernley etal. (1987) 
Reid et al. (1988) 
Whitelock (1992) 
Dwek et al. (1995) 

Binnev et al. 

Feast (1997) 
Staneketal. (1997) 
Freudenreich (1998) 



Paczvnski & Stanek (1998) 
Udalski (1998) 
Udalski (1998) 
Sevenster et al. (1999) 
Bissantz & Gerhard (2002) 



Eisenhauer et al. (2003) 
Robin et al. (2003) 
Merrifield (2004) 



8.0: 
7.1 : 

9.1 



0.65 
1.5 



8.1 ± 0.4 



8.4 ± 0.4 
8.1 ±0.15 
8.1 ± 0.06 



7.94 ± 0.42 



1:0.25:0.25 
1:0.33±0.11:0.23= 



0.08 



1.9 100 1:0.6:0.4 



1:0.43:0.29 



2.6 



2.5 

2.8 100 l:(0.3-0.4):0.3 



1:0.27:0.27 
1:0.6:0.4 



Babusiaux & Gilmore (2005) 7.7 ± 0. 15 

Eisenhauer et al. (2005) 7.62 ± 0.32 
Groeneweeen & Blommaert (2005) 8.8 ± 0.4 
Looez-Corredoira et al. (2005) 

Nishivama et al. (2006) 7.51 ± 0.10 ± 0.35 

Looez-Corredoira et al. (2007) 

Rattenburv et al. (2007) 

Ghez et al. (2008) 8.0 ± 0.4 

Gillessen et al. (2009) 8.33 ± 0.35 



1:0.5:0.4 



1:0.35:0.26 



45 

20: 

20 



10 



20-30 



44 

20- 



25 



11.1 

25 



:0.7 



22 ± 5.5 



47 

20- 

43 

24- 



35 



27 



RR Lyrae stars 

H2O maser spots 

Mira variables 

COBEfDIRBE 

surface brightness map 

COBE/DIRBE 

surface brightness map 

RR Lyrae stars 

red clump stars 

DIRBE full-sky surface 

brightness map 

red clump stars 

RR Lyrae stars 

red clump stars 

OH/IR stars 

COBEfDIRBE 

L-band map 

stars orbiting black hole 

Hipparcos data 

Hi gas and COBEfDIRBE 

surface brightness map 

red clump stars 

stars orbiting black hole 

Mira variables 

2MASS star counts 

red clump stars 

red clump stars 

red clump stars 

stars orbiting black hole 

stars orbiting black hole 



2.3. Star formation rate and metallicity distribution 

Besides the uncertainty in the literature on the parameters 
describing the GB geometry and its bar, the characteristics of 
the stars located in the GB are being debated. 

Lately, studies show traces of an intermediate age population 
in the GB. Before, it was believed that the GB consisted only of 
an old population (more than 10 Gyr). Recent studies still reveal 
that the GB is dominated by this old population, but traces of a 
smaller interme d iate age population can no longer be ignored. 
iHoltzman et ail £[993 ) studied a field in Baade's Window 
and concluded that, based on the luminosity function, there 
exist not only old stars in the GB, but also a significant, 
although unqu alified, number of intermediate age stars (less 
than 10 Gyr). lOrtolani et all (1 19951) conclude that there is no 
age difference between the majority of bulge stars and the 
halo globular clusters and no more than ~ 10% of the bulge 
popu lation can be repre sented by i ntermediate age stars (see 
also lOrtolanietal.ll200lL ZRO2003). lFeltzing & Gilmorel lfeoOO) 
found no significant young stellar population in the GB, but 
emph asise that it is s ti ll poss ible to have an age range of several 
Gyr. Ivan Loon etall (I2003L hereafter vLGO2003) found that 
in addition to the dominant old population (> 7 Gyr), there 
is also an intermediate-age population (~ 200 Myr - 7 Gyr) 
and possibly even younger (< 200 Myr) stars are found across 
the inner bulge (although the latter could be attributed to 
foreground stars). Concerning an intermediate age population, 
iGroenewegen & Blommaert (2005) found that the Mira vari- 
ables Jn_fhe_s^udi^dOGLEbulge fields have ages of a few Gyr. 
Also lUttenthaler et all d2007l) found some AGB stars that might 



be tracers of a younger population: some of their selected AGB 
stars in the GB show Technetium, which indicates that these 
AGB stars must originate from a younger population. 



The detection of a metallicity spread or absence of it and 
the metallicity distribution reveals information on the formation 
history. During the last decade, there were several papers 
studying the metallicity in the GB. iMc William & Rich! $1994) 
obtained for 14 M giants in the GB high resolution spectroscopy. 
They found a mea n metallicity of <[Fe/H]>= -0.19 ± 0.02. 
Sadl er et ail 0996) studied K and M giants in Baade's Window 
and found a mean abund ance of <[Fe/H]>= -0.11 ± 0.04. 
iFeltzing & Gilmorel (l2000h studied HST images and found that 
the metallicity of the bulge is equal to that of the old disk and 
that there is on ly a marginal e videnc e for a central metallicity 
gradient. Also Ramirez et al.] (120001) studied M giants in the 
GB. They found a mean metallicity <[Fe/H]>= -0.21 + 0.30 
and no evidence for a metallicity gradient along the minor or 
major axes of the inner bulge (R < 560 pc). vLGO2003 found a 
mean metallicity [M/H] ~ 0.5, but also stars with a metallicity 
of +0.5 and -2 are common, but not dominant. The old stars in 
the population (> 7 Gyr) tend to have higher metallicities (see 
their Fig. 18). They propose that the metallicity distribution of 
the old stars might be bimodal: one component of super-solar 
metallicity and another of subsolar metallicity. ZRO2003 found 
that most of the bulge stars have a metallicity [M/H] between 
+0.1 and -0.6 with a peak at -0.1. 
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3. TRILEGAL 

In this paper, we use the code TRILEGAL ( TRIdimensional 
modeL of thE GALaxy, Girar di et al.l Hp05), a population 
synthesis code to simulate the stellar content towards any 
direction on the sky. In this paper the model will be used to 
compute colour-magnitude diagrams (CMD) towards the GB. 
This means that we simulate the photometric properties of stars 
located towards a given direction (I, b) and complete down to a 
given limiting magnitude. 

Fig. 1 in Girardi Tt al.l d2005l) shows the general scheme of 
the TRILEGAL code. As is shown in their Fig. 1, the input to 
the model consists of four main elements: 

- a library of stell ar evolutionary track s. It is essentially 
composed of the Girar di et alj d2000l) tracks with some 
updates at low metallicities (Girardi 2002, unpublished), 
and now complemented with the TP-AGB models by 
Mari go & Girardil {2007). These latter are expected to de- 
scribe the TP-AGB phase with much m ore realistic li f etimes 
and luminosities than in the oris inal iGirardi et all ([2000) 
isochrones. 

- a library of synthetic spectra, which is used to compute 
bolometric corrections and extinction coeffic ients for the 
simula ted stars. They are described in detail by Girar di et alj 
(l2002h . 

- the instrumental setup, describing the settings of the tele- 
scope to be used, e.g. filters, detectors, and the effective sky 
area to be simulated. 

- a detailed description of the Galaxy components, be- 
ing the Galactic thin and thick disk, the Halo and the 
Bulge. For each component a Star Formation Rate (SFR), 
Age-Metallicity Relation (AMR), and space densities are 
assumed. These parameters are different for every Galaxy 
component, therefore the Galactic components are treated 
separately by the code. Also the Initial Mass Function (IMF) 
and interstellar absorption belong to this group of input 
parameters. 

The TRILEGAL code itself performs a Monte Carlo 
simulation to generate stars in accordance with the input (see 
above). The SFR, AMR and IMF define a stars age, metallicity, 
and mass. An interpolation in the grids of evolutionary tracks 
results in the absolute photometry, which is then converted 
to apparent magnitudes based on the bolometric corrections, 
distance modulus, and extinction. 

The output of the TRILEGAL code is a catalogue of stars 
that contains for each star the Galactic component in which 
they are located, stellar age, metallicity [M/H], initial mass, 
luminosity L, effective temperature r e ff, gravity g, distance 
modulus m-Mo, visual extinction Ay, core mass and the perfect 
photometric data in each of the wanted filters. To be able to 
compare the output in CMDs with observational data one has 
to perform some extra calculations which are not included in 
the TRILEGAL code such as adding photometric errors, sat- 
uration of the bright stars in real observations, and completeness. 

A more detailed description of the TRILEGAL code and the 
parameters used concerning the disk and halo can be found in 



IGirardi et al.1 (120051) . 

Extinction is calculated by the TRILEGAL code for each 
object separately based on the input Ay value and the distance 
modulus of the objec ts. Therefore, the ex tinction varies along 
the line-of-sight (see IGirardi et alj I2005L for a more detailed 
description on the way the extinction is distributed). Some 
c hanges are made, h owever, to this part of the code with respect 
to lGirardi et alj d2005l) . Not only can we now give the extinction 
value Ay as an input parameter to the code, but a spread on the 
Ay has also been used in this work. The extinction value s and 
their s igma used can be found in Tableland are based on Sumi 
(12004 . 

An important point is that we have used T RILEGAL with 
the o riginal filter cur ves provided for 2M ASS dSkrutskie et al.l 
2006) and OGLE-II dUdalski et all 1 19971) . so simulating as far 
as possible the original photometric system^]. The bolometric 
co rrections for 2MASS have already been described in detail 
bv iBonatto et all d2004l) . Those for OGLE-II have been newly 
computed, and apart from being used in TRILEGAL have 
also been incorporated into the Padova database of isochrones. 
Extensive tables of bolometric corrections and isochrones in 
the OGLE-II system are now made available at the websites 
http : //pleiadi . pd . astro . it/isoc_photsys . 82 (static 
tables) and http://stev.oapd.inaf.it/cmd (interactive 
web form). 

Furthermore, extinction coefficients in the several 
2MASS and OGLE-II passbands h a ve be en computed for 
a G2V star using the ICardelli et alJ dl989l) extinction curve 
with Ry = 3.1. The extinction coefficients turn out to be 
Auo^lAy = 1.524, A BoGL JA v = 1.308, Ay OGL JA v = 1.021, 
Ai^JAy = 0.571, A JwA JA v = 0.290, A HlMASS /A v = 0.183, 
and A Ks2MAS JA v = 0.118. 

These numbers would be suitable to simulate extinction in 
the nearby disk, but do not apply to a bulge with an anomalous 
ratio of total to selective extinction (see Udalski 2003). Since the 
reasons for such anomalous ratios are still not understood, one 
does not know really how to properly simulate extinction for the 
bulge. Our choice has been very pragmatic: we assume that the 
generalised Carde liTet al.l (119891) extinction curve still applies, 
and just change its input Ry v alue until the mean Ryi = 1.964 
value for the bulge (Sumi 2004) is recovered. Therefore, we find 
that the anomalous extinction is reproduced with a Ry - 2.4, 
which produces A UoGLE /A v = 1.715, A BoGLE /A v = 1.386, 
A VooL JA v = 1.017, Aj QGLE jAy = 0.506, Aj 1Mh JA v = 0.256, 
^massMv = 0.162, A Ks2MAS JA v = 0.101. These values are 
assumed in this work. 



4. Observational Data 

4. 1 . Two Micron All Sky Survey - 2MASS 

The Two Micron All Sky Survey (2MASS) project is a col- 
laboration between The University of Massachusetts and the 
Infrared Processing and Analysis Center (JPL/Caltech). The 

1 Our experience is that the simulation of the 2MASS photometric 
system presents zero-point offsets of just a few 0.01 mag. Zero-point 
errors of similar magnitude are just expected in synthetic photometry 
applied to present-day libraries of synthetic spectra, and are also likely 
present in our simulations of the OGLE-II system. 
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Table 2. Properties of the selected fields. 



Ogle field 


e 


b 




O-Ay" 


# stars 


(buLsc) 








OGLE 


2MASS 


2 


2.23 


-3.46 


1.51 


0.024 


803 269 


26 147 


7 


-0.14 


-5.91 


1.34 


0.025 


462 748 


19 304 


8 


10.48 


-3.78 


2.14 


0.042 


401 813 


19411 


13 


7.91 


-3.58 


2.05 


0.042 


569 850 


20 704 


14 


5.23 


2.81 


2.49 


0.048 


619 028 


24 851 


17 


5.28 


-3.45 


1.92 


0.031 


687 019 


22 241 


25 


-2.32 


-3.56 


2.33 


0.042 


622 326 


26 328 


27 


-4.92 


-3.65 


1.69 


0.025 


690 785 


23 623 


29 


-6.64 


-4.62 


1.53 


0.024 


491 941 


20 692 


40 


-2.99 


-3.14 


2.92 


0.066 


630 774 


26 169 


47 


-11.19 


-2.60 


2.58 


0.056 


300 705 


20 099 



Based on lSumil d2004l) 




resulting 2MASS project made uniformly-calibrated observa- 
tions of the entire sky (with a sky coverage of 99.998%) in the 
J (1.24 ^m), H (1-66 im\) and K s (2.16 /mi) near-infrared bands 
dSkrutskie et alj|2006l) . 

In this paper, we will use the 2MASS Point Source Catalog 
(PSC), which contains astrometry and photometry in the three 
survey bandpasses for 470 992 970 sources. For every source, 
position, magnitude, astrometric and photometric uncertainties, 
and flags, indicating the quality of the source characterisations, 
are provided. 



4.2. Optical Gravitation Lensing Experiment - OGLE 

The Optical Gravitational Lensing Experiment (OGLE) aims at 
detecting dark matter using micro lensing phenomena. Crowded 
regions are a suited place to conduct this survey. The Large 
Magellanic Cloud (LMC) was the first target in the second 
phase. Later on additional targets were added among 1 1 deg 2 in 
theGB. 



In this paper, we will use data from lUdalski et alj (j2002) 
from OGLE-II who present the VI photometric maps of the GB 
as a natural by-product of large micro lensing surveys. These 
maps contain the mean V/-photometry and astrometry of about 
30 million stars, covering about 1 1 deg 2 in different parts of the 
GB. Each field covers 14.2' x 57'. 



Fig. 1. Example of the photometric errors of the A^-band 
2MASS data for field bul.sc40. 



confidence intervals for the derived parameters. 



5.1. Method 

The output of the TRILEGAL code consists of perfect pho- 
tometric data (see Sect. [3]). In order to be able to compare 
the output of the code with real observations, we need to add 
photometric errors to the output of TRILEGAL. 
Both 2MASS and OGLE-II give individual errors on the 
observations. These individual errors are shown in Fig. Q] as 
an example for the 2MASS K s -bwd data in field bul_sc40 in 
function of the K s -bwA magnitude. Using this data a density 
profile has been generated: for each bin in K s (size 0.1 mag) 
we calculated how likely it is for a star to have an error <jk s - 
Each bin in o~k, (size 0.01 mag) is then given a number between 
and 1 representing the relative amount of stars in this bin 
with respect to all stars with a similar K s -band magnitude. For 
each star in the output file of TRILEGAL a random number has 
been generated between and 1 . The assigned error is then the 
closest to cr Ks with respect to the assigned density profile. This 
method has been used to generate photometric errors in all the 
observed magnitude bands from both 2MASS as OGLE-II for 
each field individually. 



4.3. Selected fields 

The fields selected to perform this study are shown in 
Table [2] The f ields a re chosen based on the OGLE fields from 
Udals ki et alJ (2002) and are as much as possible spread in 
longitude and latitude around the Galactic centre. The data 
from 2MASS has been downloaded according to the posi- 
tions of the OGLE fields, using Gator on the IRSA web site 
(http : //ir sa . ipac . caltech . edu/) . 



5. Analysis 

When calculating a set of models of any kind, one always needs 
an objective method to determine why one model is better 
than another. One also wants to know the error bars and/or 



With the observational data and the data calculated with 
TRILEGAL, CMDs have been constructed. For the 2MASS data 
these CMDs are in J - K s vs K s (see Fig. [2j, for the OGLE data 
the CMDs are V — I vs / (see Fig. [3}. To be able to compare a 
CMD constructed with observational data to a CMD constructed 
with modelled data, the CMDs have been divided into regions. 
For both the 2MASS data and OGLE data, we define a disk re- 
gion and a bulge region in the dereddened CMDs as follows: 

K s>0 < 12.2 and {J - K s ) < 0.25 (for 2MASS disk) (3) 

/o < 14. 1 and (V - I) < 1 .4 (for OGLE disk) (4) 

K. s ,o< 10.7 and K sfi > -10.4(7 - K s ) + 18.9 (5) 
(for 2MASS bulge) 

7 < 12.4 and 7 > -8.2(V - 7) + 24.0 (6) 
(for OGLE bulge). 
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Fig. 2. J - K s vs K s CMD for field bul_sc40. The figure on the 
left are the observations, the figure on the right is a model with 
R = 8.6 kpc, a m = 2.7 kpc, a = 97 pc, r] = 0.68, £ = 0.30, 
4> = 14°, and f = 427.3. The black dots in the model (right 
panel) are stars that originate from the disk, the dark blue dots 
originate from the halo and the red dots are bulge stars. This 
colour code will be used in all figures of this type. The light 
grey lines indicate the boxes used to compute the model with the 
observations. 



bul_sc40 




& 2 2-1 6 8 4202-168 
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Fig. 3. Same figure and colour code as Fig.|2]for the same model 
but now for a V — / vs / CMD. The additional box which is not 
seen in Fig.|2]is used to compare the red clump stars. 



For the OGLE data we defined additionally a clump region 
which contains the red clump stars, as 16.1 < Iq < 12.6 and 
0.8 < (V-/)o < 1.8. 

The regions are defined on a model with no extinction. 
Based on the assumed extinction for a field, these boxes are also 
reddened, in this way, these boxes include for every field the 
same part of the CMD. The boxes are defined in such a way 
that there is only a very small contamination of disk stars in the 
bulge region and the other way around. Therefore, if we want 
to study the bulge parameters, the only part of interest is the 
bulge region. The stars in the selected regions, are used to create 
histograms. All histograms were created in one magnitude 



Fig. 4. Example showing how the red clump stars are fitted using 
the method described bv lStanek & Garnavichl (Qj?98). The black 
histogram are the stars selected from the "red clump box" in 
the V — I vs I CMD and the grey line is the Gaussian fit with 
underlying second degree polynomial. 



(K s for 2MASS and / for OGLE) and not in the colours. The 
histograms for the disk and bulge region, which are probably 
not affected by incompleteness, are used to compare the number 
counts in each bin. 

The histogram for the stars in the red clump region are 
not used for number counts, but to derive the red clump mag- 
nitude both from observations and modelled data. To derive 
the Ip.m magnitude for the red clump, the method described in 
Stanek & Garnavich (1998) has been used. This method fits a 
2D polynomial to the underlying population, with on top of this a 
Gaussian function to determine the magnitude of the red clump. 
The formula is given by: 

«(/()) = a + b(k - 7 ,m) + c(I - /o,m) 2 

x N RC { (k-k, m ) 2 ) n , 
cr RC V27r I 2o- z RC ) 

with the first tree terms corresponding to the underlying 
population with coefficients a, b, and c for the second degree 
polynomial. The last term refers to the red clump population 
itself, with parameters corresponding to a Gaussian fit: scaling 
factor N rc , sigma cr RC and peak of the observed red clump 
stars 7o,m- Figure |4] shows a fit using this method to the /-band 
observations in the bul_sc40 field. Using this fit to the red 
clump, we will compare the peak position of the red clump stars 
in the observations to the peak position of the red clump in the 
model. The expected magnitude of the red clump Mf c is often 
compared to the peak position in the magnitude distribution of 
the red clump in order to estimate the distance to these stars. 



5.2. Poisson distributed maximum likelihood 

For the stars in the selected disk and bulge areas (see Fig. [2] 
and Fig. [3} the constructed histograms are used to compare 
number counts. The typical errors on the amount of stars in 
these histogram bins are Poisson errors (~ with «, the total 
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amount of stars in bin /). Therefore the selected test to perform 
on these histogram is a log likelihood test for Poisson statistics 
dEidelman et alj|2004l) : 

N , v 

- 2 = 2 g \Vi(0) - m + m In -^_J . (8) 

In this formula 9 is the set of unknown parameters one wants to 
derive, n = (m, n%, ... , is the data vector containing the 
observations with yV the number of bins in a histogram, v are the 
expected values, which are derived from the histograms of the 
modelled data and are therefore dependent on 9. When «, = 0, 
the last term in Eq. © is set to zero. 

For each field, there are 4 histograms (one for the disk re- 
gion and one for the bulge region, both for the comparison with 
2MASS and OGLE). When determining the parameters for the 
bulge, we only use the histograms based on the stars in the 
bulge boxes. This means that for each field, 2 histograms re- 
main. Table|2]shows that we selected 1 1 fields. Including all this 
gives us a value I for each model, which we want to minimise: 

11 2 

l = YjYj-2ln(A kj (9)). (9) 

7=1 k=\ 

5.3. Gaussian distributed log-likelihood 

In addition to a box for the disk and bulge stars, Fig.[3]also shows 
a box in which the red clump stars are expected to be. Equation|7] 
is fitted to the number counts in this box, as illustrated in Fig. |4] 
The same is done for the model stars. The parameter of interest 
is the mean magnitude of the RC, /o, m - 

When evaluating calculated models based on the position of 
the red clump, we cannot use the Poisson distributed maximum 
likelihood (see Sect. 15. 2\ . In this case the errors are not Poisson 
distributed, but normal N(0, o 2 ). To determine the appropriate- 
ness of a model concerning its positi on of the red c lump, the 
log-likelihood function I can be used dDecin et alj|2007t) : 



VI. , . , . /r— . 1 /^obsj _ ^model,j(#)\ 2 

Z = ^ln(t7-) + ln(V27r)+-l —J ■ 

7=1 ' 



To find the best model, this equation needs to be minimised. As 
the standard deviation cr we selected 0.025 mag, which is the 
half of the bin size used to construct the histograms containing 
the red clump stars (see also Sect. 16.31 ). 
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Fig. 5. The left figure shows the SFR normalised to 1 and the 
right figure shows the metallicity distribution normalised to 1 
(based on ZRO2003). 



our results as the stars fainter than this limit have no relevance 
to our analysis (see Sect. 15. U . Because these areas are chosen 
well above the completeness limits for both the 2MASS and 
OGLE data, we did not include completeness, again to speed up 
calculations. 

The model calculations and minimisation procedure have 
been performed using Python. The program has been written 
as such that the 1 1 different fields, that need to be calculated to 
evaluate one model (see Sect. 14. 3K were calculated simultane- 
ously over 11 different CPU's. We could make use of 18 Dell 
PowerEdge PC's with a total of 54 processors of 3400MHz, 
which makes it possible to evaluate about 5 models at the same 
time. In total we evaluated ~ 10 5 models in a time span of about 
a year. 



5.5. Uncertainties on the model parameters 

As explained by iDecin et aL I (120071) the log-likelihood func- 
tion can be used to estimate the uncertainties on the model 
parameters. A model with a likelihood value I is equally good 
(10) as the best model if min(Z) + N cl i t > I. These critical values 
Merit follow a y^-distribution with p the degrees of freedom 
(EidelmanetinEQQl |De 



ecin et al.1 120071). A table with t hese 



values for p e [1,20] can be found in Decin et ail (l2007l and 
references therein). 



5.4. Minimisation procedure 

The minimisation procedure used in this work is the Broyden- 
Fletcher- Goldf arb-Shann o (BFGS) method ferovdenl Il970t 
iFletcherl 119701: iGoldfarbl Il970t IShannol [1970J) which can be 
used to solve a non-linear optimisation problem. By analysing 
gradient vectors, the method constructs an approximated 
Hessian matrix allowing a quasi-Newton fitting method to move 
towards the minimum in parameter space. The parameters were 
normalised in the minimisation procedure so they have an equal 
influence on the choice of the next set of parameters. 

To significantly decrease the amount of computing time, 
the simulations made during the minimisation runs, we adopt a 
shallower limiting magnitude, at I = 18. This does not influence 



6. Results 

The SFR and AMR chosen to start the minimisation process 
with are based on ZRO2003 and are shown in Fig. [5] The left 
figure shows a bulge formed by a single star burst about 10 Gyr 
ago and the right figure indicates that the majority of the stars 
have a metallicity slightly lower than solar (Z Q = 0.019). 

The start values for t he minimisation pro cedure have been 
chosen in accordance with Binnev ^t al.l(l 19971 see also Table[TJi: 
R = 8.0 kpc, a m = 1.9 kpc, a = 100 pc, 77 = 0.60, ( = 0.40, 
(f> = 20°, and fo = 624.0. The output of this minimisation 
procedure is a model with following parameters: /?o = 8.6 kpc, 
a m = 2.7 kpc, a = 97 pc, r\ = 0.69, £ = 0.30, 4> = 20°, and 
fo = 422.3 (see also TablefJ) model 1). 
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Table 3. Overview of the different models described in the text. 



Number 


[kpc] 


[kpc] 


a 
[pc] 


1 


£ 


4> 

n 


fo 


SFR 


[Fe/H] 




'bulge 


'rc 


1 


8.6 


2.7 


97 


0.69 


0.30 


20 


422.3 


ZRO2003, 10 Gyr 


ZRO2003 




2 185 


10.6 


2 


8.6 


2.7 


97 


0.69 


0.30 


20 


422.3 


ZRO2003, 10 Gyr 


ZRO2003 


+ 0.1 dex 


1 768 


9.2 


3 


8.6 


2.7 


97 


0.69 


0.30 


20 


420.5 


ZRO2003, 10 Gyr 


ZRO2003 


+ 0.1 dex 


1761 


9.2 


4 


8.6 


2.7 


97 


0.69 


0.30 


20 


422.3 


ZRO2003, 10 Gyr 


ZRO2003 


+ 0.2 dex 


1597 


8.8 


5 


8.6 


2.6 


96 


0.68 


0.30 


20 


422.8 


ZRO2003, 10 Gyr 


ZRO2003 


+ 0.2 dex 


1576 


8.4 


6 


8.6 


2.7 


97 


0.69 


0.30 


20 


422.3 


ZRO2003, 10 Gyr 


ZRO2003 


+ 0.3 dex 


1519 


9.2 


7 


8.6 


2.7 


97 


0.68 


0.30 


14 


427.3 


ZRO2003, 10 Gyr 


ZRO2003 


+ 0.3 dex 


1442 


13.3 


8 


8.6 


2.7 


97 


0.69 


0.30 


20 


422.3 


ZRO2003, 10 Gyr 


ZRO2003 


+ 0.4 dex 


1597 


8.7 


9 


8.6 


2.8 


97 


0.68 


0.28 


20 


421.6 


ZRO2003, 10 Gyr 


ZRO2003 


+ 0.4 dex 


1 560 


8.3 


10 


8.6 


2.7 


97 


0.69 


0.30 


20 


422.3 


ZRO2003, 10 Gyr 


ZRO2003 


+ 0.5 dex 


1768 


8.6 


11 


8.6 


2.7 


97 


0.69 


0.30 


21 


425.5 


ZRO2003, 10 Gyr 


ZRO2003 


+ 0.5 dex 


1717 


8.5 


12 


8.6 


2.7 


97 


0.68 


0.30 


14 


427.3 


ZRO2003, 7 Gyr 


ZRO2003 


+ 0.3 dex 


3517 


18.0 


13 


8.7 


2.6 


96 


0.68 


0.28 


14 


419.0 


ZRO2003, 7 Gyr 


ZRO2003 


+ 0.3 dex 


1516 


20.4 


14 


8.6 


2.7 


97 


0.68 


0.30 


14 


427.3 


ZRO2003, 8 Gyr 


ZRO2003 


+ 0.3 dex 


2247 


17.0 


15 


8.7 


2.6 


95 


0.68 


0.31 


15 


406.0 


ZRO2003, 8 Gyr 


ZRO2003 


+ 0.3 dex 


1432 


9.7 


16 


8.6 


2.7 


97 


0.68 


0.30 


14 


427.3 


ZRO2003, 9 Gyr 


ZRO2003 


+ 0.3 dex 


1735 


13.0 


17 


8.6 


2.6 


96 


0.69 


0.29 


14 


427.5 


ZRO2003, 9 Gyr 


ZRO2003 


+ 0.3 dex 


1433 


15.1 


18 


8.6 


2.7 


97 


0.68 


0.30 


14 


427.3 


ZRO2003, 1 1 Gyr 


ZRO2003 


+ 0.3 dex 


1657 


7.1 


19 


8.6 


2.7 


97 


0.68 


0.30 


15 


434.2 


ZRO2003, 11 Gyr 


ZRO2003 


+ 0.3 dex 


1576 


10.8 


20 


8.6 


2.7 


97 


0.68 


0.30 


14 


427.3 


ZRO2003, 12 Gyr 


ZRO2003 


+ 0.3 dex 


1851 


25.1 


21 


8.6 


2.7 


97 


0.68 


0.30 


13 


427.5 


ZRO2003, 12 Gyr 


ZRO2003 


+ 0.3 dex 


1 564 


11.9 


22 


8.6 


2.7 


97 


0.68 


0.30 


14 


427.3 


ZRO2003, 13 Gyr 


ZRO2003 


+ 0.3 dex 


2 342 


7.4 


23 


8.5 


2.8 


97 


0.68 


0.31 


16 


432.3 


ZRO2003, 13 Gyr 


ZRO2003 


+ 0.3 dex 


1 696 


11.5 


24 


8.7 


2.6 


95 


0.68 


0.31 


15 


406.0 


VLGO2003 


ZRO2003 


+ 0.3 dex 


2 190 


18.7 


25 


10.7 


3.1 


86 


0.73 


0.30 


8 


455.0 


vLGO2003 


ZRO2003 


+ 0.3 dex 


1747 


120.9 


26 


8.7 


2.6 


95 


0.68 


0.31 


15 


406.0 


vLGO2003, 


ZRO2003 


+ 0.3 dex 


2112 


18.9 


















older than 0. 1 Gyr 










27 


9.3 


2.8 


92 


0.66 


0.30 


17 


394.0 


VLGO2003, 


ZRO2003 


+ 0.3 dex 


1466 


11.4 


















older than 0. 1 Gyr 










28 


8.7 


2.6 


95 


0.68 


0.31 


15 


406.0 


VLGO2003, 


ZRO2003 


+ 0.3 dex 


3190 


12.0 


















older than 1.0 Gyr 










29 


9.3 


2.8 


92 


0.66 


0.30 


17 


394.0 


VLGO2003, 


ZRO2003 


+ 0.3 dex 


2749 


17.4 


















older than 1.0 Gyr 










30 


8.7 


2.6 


95 


0.68 


0.31 


27 


406.0 


ZRO2003, 8 Gyr 


ZRO2003 


+ 0.3 dex 


1426 


7.3 



The constructed histograms for this model are shown for 
field bul_sc40 in Fig. [6] As mentioned in Sect. 15.11 we only 
used the histograms constructed with the / and K s magnitude 
to compare the model with the observations. The histograms 
in V — I and J - K s illustrate the fit between the model and 
the data for the colours. The minimisation procedure uses only 
the stars in the bulge box. The stars in the disk box were used 
prior to the minimisation process to set the scale parameter for 
the disk stars and to check the contamination of disk and halo 
objects in the "bulge box". For field bul_sc40 adopting the same 
model as in Fig. [6] we find that out of a total of 2731 objects 
in the bulge box for the 2MASS data only 127 originate from 
the disk and halo. Concerning the 2110 objects in the bulge box 
for the OGLE data, only 105 originate from the disk and halo. 
Since this is already a very modest fraction (about 5% for both 
CMDs), possible errors in our representation of the disk and 
halo geometry are expected to have a negligible impact on the 
star counts inside the bulge box. 

Fig.|6]shows the distribution of stars for the observations and 
the model for the selected bulge stars for field bul_sc40. Also 
for all the other fields this model overestimates the amount of 



bulge stars in the 2MASS data and underestimates the amount 
of bulge stars in the OGLE data, therefore this model can still 
be improved, but not by adapting the parameters already in the 
minimisation procedure. 



6. 1 . Varying the metallicity distribution 

One of the parameters not in the minimisation procedure is the 
metallicity. The metallicity distribution in ZRO2003, which 
has been used for the previous model calculations, is based on 
photometric results and is slightly less metal-rich in comparison 
with meta llicity distributions derived f rom spectros c opic r esults 
from e.g. iMc William & Richl (1 1994 . ISadler et alJ (l996), and 
iRamfrez et al.1 d2000l) (see Fig. 14 in ZRO2003). Therefore, we 
shifted the metallicity distribution towards a more metal-rich 
population. Already when we shift [Fe/H] with 0.1 dex our fit 
improves (see model 2 in Table [3]). The metallicity distribution 
has been shifted with +0.1, +0.2, +0.3, +0.4, and +0.5 dex 
(see respectively models 2, 4, 6, 8, and 10 in Table [3j. We 
obtained the best fit with a shift of +0.3 dex for [Fe/H]. New 
minimisation procedures have been set up using these shifted 
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Fig. 6. Histograms showing the distribution of selected bulge 
stars for field bul_sc40. The black line represents the observa- 
tions, the red line the model. The left panel shows the 2MASS 
data, the right panel the OGLE data. The model has the fol- 
lowing parameters: Ro = 8.6 kpc, a m = 2.7 kpc, ao = 97 pc, 
r, = 0.69, ( = 0.30, 4> = 20°, and f = 422.3. The star formation 
rate is a 10 Gyr burst and the metallicity distribution is based on 
ZRO2003. 
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Fig. 7. Histograms showing the distribution of selected bulge 
stars for field bul_sc40. The black line represents the observa- 
tions, the red line the model. The upper panels shows the 2MASS 
data, the lower panels the OGLE data. The model has the fol- 
lowing parameters: Rq = 8.6 kpc, a m = 2.7 kpc, ao = 97 pc, 
r, = 0.68, C = 0.30, 4> = 14°, and f = 4273. The star formation 
rate is a 10 Gyr burst and the metallicity distribution is based 
on ZRO2003 but shifted with 0.3 dex. (This figure is also avail- 
able in electronic form, showing histograms for all the modelled 
fields.) 



metallicity distributions in order to test the effect of the change 
in metallicity distribution on the derived parameters. The 
minimisation procedure was started from the result of the 
previous minimisation process (see model 1 in Table [3}. The 
results are listed in Table [3] (see models 3, 5, 7, 9, and 11). 
Reminimising the models revealed no differences concerning 
the metallicity distribution shift, again the best results are 



found using a shifted metallicity distribution of +0.3 dex with 
respect to ZRO2003. Concerning the parameters describing the 
geometry of the bulge and its bar, only and /o change. These 
changes are small and fall within the expected error bars (see 
Sect. |6J]l- Fig. 14 in ZRO2003 shows that a shift with 0.3 dex 
in [Fe/H] is consistent with the noticed differences between the 
metallicity distribution they derived and me tallicity distributions 
derived based on spectroscopy. In fact, IZoccali et al.l ((2008) 
have recently rederived the bulge metallicity distribution, based 
on high resolution spectroscopy for a large sample. They find 
metallicities systematically higher than in ZRO03 by about 0.2 
dex, which is perfectly in line with what we obtain in this paper. 



6.2. Varying the age distribution 

The star formation rate is also a input to the model (see 
Sect. I2.3l l. Therefore we also tested other SFR's. First we start 
with shifting the star burst to a different age in steps of 1 Gyr 
from a star burst 7 Gyr ago to a star burst 13 Gyr ago. Table [3] 
gives an overview of these models when no new minimisation 
procedure is started (see models 12, 14, 16, 18, 20, and 22). 
Amongst these models, there is no better model than model 7. 
When a minimisation procedure is set up, using the different 
age distributions and the parameters describing the geometry 
as derived in model 7 as start values, we find models 13, 15, 
17, 19, 21, and 23 (see Table [3]). Three of these models are 
equally good based on /bulge: model 15 with a star burst of 8 Gyr 
is the so-called best model, but models 17 (star burst of 9 Gyr 
ago) and 7 (star burst 10 Gyr ago) are equally good and one 
cannot distinguish between these three models. Concerning 
the parameters describing the geometry, the differences seen 
between these three models, fall within the expected error bars 
(see Sect. 16.3b . 

Changing the age of the star burst is not the only possibility 
to alter the SFR. Sect. 12.31 describes that most likely there are 
also intermediate age stars in the GB and maybe also some 
younger stars. In order to check these possibilities, we used the 
SFR given in vLGO2003 (see their Fig. 28). The metallicity dis- 
tribution used is the metallicity distribution of ZRO2003 shifted 
with 0.3 dex. Young stars as well as intermediate age stars are 
included in the SFR by vLGO2003. To check the possibility 
of an intermediate age population on top of the dominant old 
population, the age distribution of vLGO2003 has been used, 
but without the stars younger than 0. 1 Gyr and a SFR without 
the stars younger than 1 Gyr. The metallicity distribution used 
in these scenarios is again the shifted distribution found by 
ZRO2003. The results are again shown in Table [3] 

Model 24 in Table[3]shows that we do not get a significantly 
better fit if we use the SFR by vLGO2003 and the shifted 
metallicity distribution by ZRO2003. Because the SFR has 
changed drastic, and therefore also the colour distribution of the 
stars in the resulting models, a new minimisation procedure is 
necessary to check if no better model can be found using the 
SFR by vLGO2003. The previous best model found (model 15), 
has been chosen to start the minimisation procedure with. This 
resulted in the following model: Ro = 10.7 kpc, a m = 3.1 kpc, 
a = 86 pc, j] = 0.73, £ = 0.30, cf> = 8°, and / = 455.0 (model 
25 in Table [3}. All parameters except Rq and fo fall within the 
expected error bars of the best model with a star burst of 10 Gyr 
as SFR (see Sect. 16.3b . The differences seen between Ro and fo 
are huge and do not fall within the previously defined error bars. 
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'bulge Additionally, shows that this model is not significantly 
better than model 15. 

Also for the two SFR's including only intermediate age 
stars based on the SFR by vLGO2003 a new minimisation 
procedure was set up. Using the SFR by vLGO2003 for stars 
older than 0. 1 Gyr, the minimisation procedure ended with the 
following model: Rq = 9.3 kpc, a m = 2.8 kpc, «o = 92 pc, 
r] = 0.66, £ = 0.30, 4> = 17°, and f = 394.0. All parameters 
found using this SFR are within the error bars derived on 
model 15 (see Sect. 16.31 ). Using the SFR for stars older than 
1 Gyr, the minimisation procedure ended with the same model: 
Rq = 9.3 kpc, a m = 2.8 kpc, a = 92 pc, 77 = 0.66, f = 0.30, 
(p = 17°, and fa = 394.0, which then of course falls also within 
the error bars of model 15. Based on /bulge, these two later 
models are not as good as model 15. 



6.3. Error bars 

To derive the error bars on these model parameters, we used the 
method described in Sect. 15.51 For each parameter we explored 
the parameter space around the value found by the minimisation 
procedure and rounded the found parameters within the found 
error bars. All models that fall within the critical value of 
N c „t = 14.07 from our model are equally good. This results 
in the following error bars for model 7: Rq = 8.6±2}q kpc, 
a m = 2.7+0 04 kpc ^ a(j = 97± i.i pC; ^ = 0.68+0 01, ( = o.30+o o| , 

0=14°±« and/ = 427.3±^. 

These error bars have been used in the previous section to 
determine whether parameters between different models fall 
within the same range. As model 15 was found to be the best 
model, we also calculated the error bars for this model. We 
obtain: Rq = 8.7±°,42 kpc, a m = 2.5±J' ™ kpc, «o = 95±j 2 5 g pc, 
r, = 0.68+0;04, f = 6.31±j$J, cf> = 15±%*. and/ = 406.0±^. 
The error bars on this model are somewhat larger than the error 
bars found for model 7. 

Based on the outcome of the minimisation procedures, the 
error bars on the determination of the magnitude of the red 
clump stars in the /-band have been determined. This resulted 
in an average error of 0.03 mag which is consistent with the cho- 
sen value in Eq. ( fTOb . 

7. Red Clump stars 

The red clump stars can be used as a distance indicator as well as 
an estimator to the angle <p between the bar and the Sun-centre 
line (see also Sect. 12.21 1. Using the red clump stars, a new 
minimisation procedure as described in Sect. 15 .41 was set up. In 
this scenario, there are only two variable parameters: Rq and <p, 
the other parameters have been kept constant (a m = 2.6 kpc, 
a = 95 pc, 77 = 0.68, £ = 0.31, and f = 406.0). An age of 
8 Gyr and a metallicity distribution based on ZRO2003 and 
shifted with 0.3 dex has been used (which is consistent with the 
best model found, 15). The minimisation process started with 
Rq = 8.7 kpc and </> = 14° and ended with Rq = 8.7 pc and 
<p - 26.7°. This model is also listed in Table|3] 

Using the method to derive error bars on these parameters, 
with A^rit = 5.99, the following error bars were calculated: 
Rq = 8.7+og kpc and = 26.7°±* 5 r Therefore the model 
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Fig. 10. Peak positions Iq of the red clump stars corrected for 
extinction in function of the longitude I. 



derived based on the red clump stars for Rq and cf> falls within 
the previously defined error bars for the model derived on the 
stars located in the bulge box. Although this model has a smaller 
/bulge and /rc, it is not significantly better than the previous 
found best model (see model 15 in Table[3]) 

In Sect. [6] other SFR's have been tested (see Table [3]). Also 
for these other models, we compared the red clump stars. Based 
on the red clump stars, Table[3]shows that the models with a star 
burst of 10 Gyr as SFR but with shifted metallicity distributions 
are equally good, except for model 7 which has a significantly 
worse /rc - There is based on the critical value N a -\ t no difference 
between these models. Concerning the models with a SFR based 
on vLGO2003 the comparison with the red clump stars shows 
that our fit could not be improved adapting this SFR. 

The results of this minimisation, and in particular the Rq 
value, are clearly associated to the set of stellar models that 
we used. Models with a different absolute magnitude of the 
red clump, for insta nce, would lead to a d i fferen t GC distance. 
iGirardi et all dl998l) and lGirardi & Salarisl d200ll) claim the set 
of stellar models we are using dGirardi et al.ll2.000l) do reproduce 
the very accurate Mf c value derived from Hipparcos; however 
their claims are tight to the specific choi ce of SFR (from 
iRocha-Pinto et al.1 lioOObl) and AMR (from iRocha-Pinto et al.1 
l2000al) employed in modelling the Solar Neighbourhood. Future 
revisions of these functions may lead to different results. 

Several other works have already tried to access the bulge 
orientation and GC distance using red clump stars, employing 
different data, extinction corrections, and assumptions about the 
intrinsic magnitude difference between the local (Hipparcos) 
and the bulge red clump (i.e. the population correction, 
AM RC = M RC (Hipp) - M RC (galaxy)). ^From the entries in 
Table [1] one notices that our best Rq va lue agrees with the 
distan ce d eterminations to the bulge f rom iPaczvnski & Stanekl 
(1 19981) and lStanek & Garnavichl (fl998). which were also derived 
from /-band data from OGL E. However, these are signific antly 
larger than those derive d from Babusia ux & Gilmord d2005l) and 
Nishivam a et al.l d2006l) using near-infrared data. 
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In order to allow any future comparison between our Rq 
determination and those from other authors, we provide below a 
representative value of the absolute red clump magnitude of our 
models. 

Since TRILEGAL simulates the entire population along a 
given line-of-sight, our simulated LFs for the red clump box 
include effects like the dispersion in the distances of bulge 
stars, geometrical factors, and the contamination by disk and 
halo stars. As a consequence, our effective Mf c values do vary 
slightly from region to region even though the bulge stellar 
population which is being simulated is everywhere the same. 
Moreover, the Galactic centre line-of-sight has not been simu- 
lated, and a representative value of Mf c has to be inferred from 
the values distributed across the bulge. Figure [10] shows the ob- 
served peak positions of the red clump, 7n,m as a function of lon- 
gitude I. These peak positions are corrected for extinction using 
the extinction values given in Table [2] The full line is a linear 
least squares fit through these peak positions and is given by 
/o, m = (14.502 + 0.018) + (-0.029 + 0.03X. Therefore the peak 
position I , m in the GC {( = 0) is given by 14.502 + 0.018. This 
is related to the absolute red clump magnitude by 



M 



RC 



Io,m-(5 log/?,,- 5). 



(11) 



Using the distance to the GC Rq = 8.7 ±042 kpc derived 
from our minimisation procedures for an age of 8 Gyr, and 
the peak position Io jn inferred for the GC, this equation 
gives Mf c = -0.196+g Jyg. This value is consistent with the 
Mf c = -0.185 + 0.016 used by iPaczvnski & Stanekl d 19981) 
and w ith the Mf c = -0.23 + 0.03 used bv lStanek & Garnavichl 
( 1998), in their GC distance determination^. The value is also 
co nsistent with the recen t determination of Mf c = -0.22 + 0.03 
by IGroenewegen I (|2008) from an analysis based on the revised 
Hipparcos parallaxes. This explains why their distance deter- 
minations are consistent with ours despite the quite different 
methods employed. 

On the other hand, our prefe rred Rq value disagre e s with 
the iBabusiaux & Gilmord d2005l) and iNishivama etall d2006l) 
distance determinations based on near-infrared data for the red 
clump. These works find Ro values close to 7.5 kpc, which means 
a disagreement at the level of 0.3 mag in distance modulus. The 
origin of this discrepancy is not understood at the moment. 

It is worth remarking that our favoured Rq value of 8.7 kpc 
comes from both the red clump and the bulge box, the latter 
consisting mainly of upper-RGB and TP-AGB stars. These stel- 
lar groups are quite independent in their photometric properties. 
Our Ro value is also consistent with the determinations based 
on the RR Lyrae and Mira variables. Although it disagrees with 
with the BH g eometrical distance derived from Eisenhaue r et alj 
d2003l 1200 5) it does agree with the recent determination by 
iGhez etalJd2008h . 



8. Summary and conclusions 

In this paper we have studied the characteristics of the Galactic 
bulge based on CMD comparisons between OGLE and 2MASS 



2 We recall that their Mf c values were derived from the Hipparcos 
sample of red clump stars with good parallax measurements, without 
applying any population correcti on between the Hipparcos and bulge 
samples. Girar di & Salarisl d200lL their Table 4) estimated that this cor- 
rection is smaller than 0. 1 mag. 



data an the resul ts from the galactic model TRILEGAL 
dGirardi et al .112005 ). Several star formation rates and metallicity 
distributions have been tested. Over all these different SFR's and 
metallicity distributions the model parameters give within the 
error bars the same results except for the angle <p between the 
Sun-centre line and the Bar. Concerning the other parameters 
we found a distance to the Galactic centre Ro = 8. 7 ijj^Z kpc, a 
scale length for the Bulge a„, of 2.5±H? kpc, an inner truncation 
length for the Bulge qq of 95+j 2 5 8 pc, and a scaling factor 
fo = 406.0+jg 7 4 r Concerning the characteristics of the Bar, we 
found the ratio of the major axis to be 1 : 0.68 : 0.3 1± '^ 
and an angle <p of 15°+[H. The largest scatter in the found 
model parameters is found for the angle (p. Nevertheless 
the found values are all rather small and around 15 ~ 20°. 
The present best model of the bulge will be provided as 
the default in the interactive web interface to TRILEGAL 
(http : //stev . oapd . inaf . it/tril egal). 

If we compare our results to the results found in the litera- 
ture and listed in Table Q] our parameters fall within the listed 
range of parameters except for the axial ratio rj. Concerning 
the distance to the GC, our v alue resembles the values found 
by Paczvnski & Stanek (1998) based on red clump stars and 
IGroenewe gen & Blommaert (200 51) based on Mira variables. 
Our results disagree with e.g. Udalski (1998) based on red 
clump stars in the GB. If we compute the mean observed m agni- 
tude /(> fro m the distance modulus liste d in lUdalskH {T998) using 
Mf c from lStanek & Garnavichl dT998) as they did, than this is 
consistent with the observed mean magnitude we derived in this 
work. The other two studies deriving the distance to the GC 
base d on red clump stars which a re in consistent with th i s wor k 
are IBabusiaux & Gilmord d2005l) and INishivama et alj J2006). 
These two studies are consistent with each other and use a very 
similar method: both define an "extinction-free" magnitude. It 
is not exactly clear why the results of these methods are not 
consistent with our results. As these studies are performed using 
/T-band observations, w e can no t comp are the expected magni- 
tude as we could do for lUdalskH (119981). For the scale leng th of 
the bulge, we fo und the same va l ue as iFreuderireichld 19981) and 
simila r values as lSevenster et al.l (1 19991) andlBissantz & Gerhard! 
(2002). The inner truncation length f or the Bulge i s simil ar to 
the two values listed in T able Q] by Binn ev et al.l d 19971) and 
iBissantz & Gerhard] d2002l) . The value 77 we found is higher 
than the values found in the literature. Concerning the other 
ratio o f the major axis { w e foun d similar resu l ts as Whitelockl 
( 1 19921) . iDwek et all d 1995b . and IStanek et alj di997|). For the 
last parameter 0, we found values consistent with the group of 
studies finding lowe r values for (see e . g. IDwek et alj 1 1995; 
i Binney et ail I l997t IStanek et all I l997t IBissantz & Gerhard ! 
2002[lRobin et al.ll2003HMerrineldll2004l: iBabusiaux & Gi lmore 
|2005[ lLopez-Corredoira et alj 120051: iRattenburv et al.l 120071) 
As already mentioned by IGroenewegen & Blommaert! P005) 
the different results found in the literature with respect to the 
angle between the ba r and the Sun-centre line (see 



e e.g. our 

results in contrast to IWhitelockl Il992t ISevenster et al.l Il999t 



I Groenewegen & Blommaert] 120051: lLopez-Corredoira et al.l 
2007) could originate from a different spatial distribution of the 
stellar population these studies trace. 

We also tested different star formation rates and metallicity 
distributions. A model with a star burst of 8 Gyr based on 
ZRO2003 gave the best results together with a metallicity 
distribution also based on ZRO2003 but shifted by 0.3 dex, 
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although we could find no clear discrepancy between this 
model and the models with a star burst 10 Gyr ago could be 
found. We could not improve our fits including intermediate 
and/or young stars if we compared both the stars we selected 
in the defined bulge box and the red clump stars. Therefore 
we conclude that based on 2MASS and OGLE data compar- 
isons with the TRILEGAL model we find best results for a 
population of 8 Gyr. Nevertheless it is still possible that there 
are also intermediate age stars located in the GB, although 
we could not trace them using our described method (se e e.g. 
iGroenewegen & Blo mmaert 2005; Uttentha ler et a 020071) . 

Using the peak positions of the red clump stars and the 
distance to the Galactic centre retrieved by our minimisation 
procedures, we found that Mf c - -0.170+^^. This value is 
consistent with recent values foun d in the literature, althou gh it 
resembles best the values found bv lGirardi & Sa laris ( 200 j]) us- 
ing red clump stars in the solar neighbourhood {Mf c = -0.171) 
or using red clump stars in Baade's Window if a-enhancement 
is taken into account (Mf c = -0.161). 
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Fig. 9. Histograms showing the distribution of selected bulge stars for all the modelled fields. The black line is the histogram for 
the observational data, the red line represents the model. For each field there are four panels, the first two panels shows the 2MASS 
data (one for star counts in J - K s and one for the star counts in K s , The other two panels show the star counts for the OGLE data 
(V — I and /). The model has the following parameters: Rq = 8.6 kpc, a m = 2.7 kpc, ao = 97 pc, r\ = 0.68, ( = 0.30, <p = 14°, and 
/o = 427.3. The star formation rate is a 10 Gyr burst and the metallicity distribution is based on ZRO2003 but shifted with 0.3 dex. 
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Fig. 9. Continued. 



